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■ Abstract. We have carried out a detailed study of the single-degenerate channel for the pro- 
genitors of type la supernovae (SNe la). In the model, a carbon-oxygen white dwarf (CO WD) 
accretes hydrogen-rich material from an unevolved or a slightly evolved non-degenerate com- 

' panion to increase its mass to Chandrasekhar mass limit. Incorporating the prescription of 
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I Hachisu et al. (1999a) for the accretion efficiency into Eggleton's stellar evolution code and as- 
■ suming that the prescription is valid for all metallicities, we performed binary stellar evolution 
' calculations for more than 25,000 close WD binary systems with various metallicities. The ini- 

tial parameter spaces for SNe la are presented in an orbital period-secondary mass (logPi,M2) 

plane for each Z. 

(-H I Adopting the results above, we studied the birth rate of SNe la for various Z via binary 

rS |i population synthesis. From the study, we see that for a high Z, SNe la occur systemically earlier 
I and the peak value of the birth rate is larger if a single starburst is assumed. The Galactic birth 



rate from the channel is lower than (but comparable to) that inferred from observations. 

We also showed the distributions of the parameters of the binary systems at the moment 



of supernova explosion and the distributions of the properties of companions after supernova 

explosion. The former provides physics input to simulate the interaction between supernova 

ejecta and its companion, and the latter is helpful to search for the companions in supernova 

" . ' ' remnants. 
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1. Introduction 

■ Although type la supernovae (SNe la) appear to be good cosmological distance indi- 
cators and have been applied successfully in determining cosmological parameters (e.g. 
n and A: iRiess et al. 19981 IPerlmutter et al. 1999) . the exact nature of SNe la is still un- 
clear, especially the progenitor model of SNe la (see the reviews by |Hillebrandt fc Niemeyer 2000[ 
Leibundgut 2000). At present, the single-degenerate Chandrasekhar model f Whelan fc Iben 19731 



Nomoto, Thielemann fc Yokoi 1984) is the most widely accepted model. In the model, 



a CO WD accretes hydrogen-rich material from its companion until its mass reaches a 
mass of ~ 1.378M0 (close to Chandrasekhar mass, [Nomoto, Thielemann fc Yokoi 1984 ), 



and then explodes as a SN la. The companion is probably a main sequence star or a 
slightly evolved star (WD-fMS). The discovery of the potential companion of Tycho's 
supernova also verified the reliability of the WD -I- MS model ( |Ruiz-Lapuente et al. 2004| 
llhara et al. 2007P . The purpose of this paper is to study the progenitor model compre- 
hensively and systematically. 
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Figure 1. The evolution of the birth rate of 
SNe la for a single starburst of 10^^ Mq for 
different metallicities with qce = 1.0. 
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Figure 2. The evolution of the birth rate 
of SNe la for a constant star formation rate 
(Z=0.02, SFR=5A/0yr-^) with acE = 1.0. 



2. Binary evolution calculation 

We use the stellar evolution code of Eggleton fl971 [TO7^|1973j] to calculate the binary 
evolutions of WD+MS systems. The code has been updated with the latest input physics 
over the last three decades ( |Han, Podsiadlowski fc Eggleton 1994|[PoTs et al. lOQSllPols et al. 1998|) . 



Roche lobe overflow (RLOF) is treated within the code described by Han et al. (2000) 



Ten metallicities are chosen here (i.e. Z =0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 



0.04, 0.05 and 0.06). The opacity tables for these metallicties are compiled by Chen & Tout (2007) 



from Iglesias & Rogers (1996) and Alexander & Ferguson (1994) 



Instead of solving stellar structure equations of a WD, we adopt the prescription in 



Hachisu et al. (1999a) on the accretion of the hydrogen-rich material from its companion 
onto the WD based on an assumption of optically thick wind (|Hachisu et al. 1996P (see 
|Meng, Chen &: Han. 2008| for the prescription in details). 

We calculated more than 25,000 WD+MS binary systems with various metallicities, 
and obtained a large, dense model grid. We summarize the final outcomes of all the binary 
evolution calculations in the initial orbital period-secondary mass (logP', Mg) planes (see 



Figs 1 to 10 in Meng, Chen & Han. 2008). The contours for SNe la are shown by solid 
lines in these figures. From these figures, we see that the contour moves from the left lower 
region to the right upper region with metallicity Z in the (logP', Afj) plane. This means 
that the progenitor systems for SNe la have a more massive companion and a longer 
orbital period for a high Z. We wrote our results into a FORTRAN code and the code can 
be downloaded on X. Meng's personal web site http://www.ynao.ac.cn/"" bps/download 
/xiangcunmeng.htm. 



3. The results of binary population synthesis 

To investigate the birth rate of SNe la, we followed the evolution of 10^ binaries for 
various Z with Hurley's rapid binary evolution code ( [Hurley et al. 2000ll2002p . The re- 
sults of grid calculations in section [2] are incorporated into the code. Since the code is 
valid just for Z ^ 0.03, only seven metallicities (i.e. Z = 0.03, 0.02, 0.01, 0.004, 0.001, 
0.0003 and 0.0001) are examined here. The primordial binary samples are generated in a 
Monte Carlo way and a circular orbit is assumed (see |Meng, Chen fc Han. 2008] for de- 
tails). The results are shown in Figs. [1] and [D From Fig.[Tl we see that most supernovae 
occur between 0.2 Gyr and 2 Gyr after star formation, and a high metallicity leads to a 
systematically earlier explosion time. The peak value of birth rate increases with metal- 
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Figure 3. The distribution of tlie radii and 
the ratios of separations to radii of compan- 
ions at the moment of supernova explosion 
for the case of Z — 0.02 and qqe ~ 1-0. 
Cross represents the potential candidate of 
the companion of Tycho's supernova, Tycho 
G dRuiz-Lapuente et al. 2004] Branch 2004)l . 
The length of the bars of the cross represents 
observational error. Solar symbol and penta- 
cle represent the main -sequence and subgiant 
companion model in [Marietta et al. (2000) [ 
respectively. 
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Figure 4. The distribution of the masses and 
the space velocities of companions in SN la 
remnants for Z = 0.02 and qce ~ 1.0. The 
cross represents the potential candidate of 
the companion of Tycho's supernova, Tycho 
G dRuiz-Lapuente et al. 2004| Branch 2004)l . 
and the length of the bars of the cross rep- 
resents observational errors. 



licity Z. However, Fig [2] shows that the WD + MS channel can only account for about 
1/3 of the Galactic SNe la observed(3-4xl0^^yr^\ |van den Bergh fc Tammann 1991 



[Cappellaro fc Turatto 1997| . Therefore, there may be other channels or mechanisms con- 
tributing to SNe la. 

We also showed the distributions of the parameters of the binary systems at the mo- 
ment of supernova explosion (Fig. [3]) and the properties of companions after supernova 
explosion (Fig. |4]) . The former may provide physics input when one simulates the inter- 
action between supernova ejecta and its companion, and the latter may help to search 
for the companions in supernova remnants. From Fig.lH we see that the properties of the 
potential companion of Tycho's supernova, Tycho G, are well consistent with our binary 
population synthesis results. We also noticed from Fig. [3] that the solar model used by 



Marietta et al. (2000) when they simulated the interaction between explosion ejecta and 



the companion is not a typical one (see the solar symbol in the figure). 

A SN la remnant is not spherically symmetric due to the companion, and then the 
spectrum of the SN la may be polarized ([Leonard &: Filippenko 2005| . Based on the 
numerical simulation of Marietta et al. (2000) and Kasen et al. (2004) we found that at 
least 75% of all SNe la may be detected by spectropolarimetry. At present, almost all 
SNe la, which are observed by spectropolarimetry, have various degrees of polarization 
signal (jLeonard et al 2005i) . 

Since there exists mass-loss from a progenitor binary system by an optically thick 
wind, the material lost from the system may become circumstellar materials (GSM). The 
GSM may be the origin of the color excess of SNe la. Based on the single degenerate 
model with optically thick wind and using a simple analytic method, we reproduced the 
distribution of the color excess of SNe la obtained from observation if a wind velocity 
of 10 km/s is adopted. If the wind velocity is larger than 100 km/s, the reproduction is 
bad. 
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4. Conclusion 

We systematically studied the single degenerate channel of SNe la, (i.e. WD+MS chan- 
nel), and showed the parameter space leading to SNe la from this channel. Our work may 
provide help for the future studies of SNe la. We also provided the properties of com- 
panions after supernova explosion, which may help to search for the companion in SN la 
remnants. We found that at least 75% of all SNe la may be detected by spectropolarime- 
try, and then in the future, statistical studies about the spectropolarimetric observations 
of SNe la may verify the reliability of the WD -|- MS channel. 
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